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Abstract 

We study an influence of nonlinear dissipation and external perturbations 
onto transition scenarious to chaos in Lorenz-Haken system. It will be show 
that varying in external potential parameters values leads to parameters domain 
formation of chaos realization. In the modified Lorenz-Haken system transitions 
from regular to chaotic dynamics can be of Ruelle-Takens scenario, Feigenbaum 
scenario, or through intermittency. 

1 Introduction 

One of the most actual tasks in the theory of nonlinear dynamical systems 
is a setting of chaotic regime generation conditions and defining possibilities 
of chaos control [l]u]- It is well known, that in physical applications for a 
transition from chaotic to periodic mode in multicomponent systems different 
mechanisms are used. For example, in laser physics they are negative feedback 
p], angle between two crystals which are entered to the Fabry-Perot cavity 
4 . full feedback intensity |1 , intensity of activating |5 (see |6 and citations 
therein). In this connection, from the theoretical point of view, an actual task 
is to establish of chaos control and to determine transition characters between 
chaotic and regular dynamics. 

The main goal of this work is to study an influence of two additional non- 
linearityes that arise up in the chaotic system as a result of different physi- 
cal processes onto transitions character between regular and chaotic dynam- 
ics. We shall consider a modified Lorenz-Haken model which self-consistently 
can describes, for example, optical bistable systems (7], systems of defects in 
solids [8j[9], etc. Due to condition of commensurability for all three modes 
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relaxation times we shall set domains of system parameters of chaos realiza- 
tion with a help of maximal Lyapunov exponents approach. We shall obtain 
two different strange chaotic attractors and set possible transitions to chaotic 
dynamics. 

The paper is organized in the following manner. In Section 2 we present a 
model of our system incorporating a nonlinear dissipation external perturbation 
terms. Section 3 is devoted to the consideration of conditions for transition to 
chaotic regime and to determination the main characteristics of strange chaotic 
attractor. The main results and prospects for the future are presented in the 
Conclusions. 



2 A model of a chaotic system 

A Lorenz-Haken model can be written in a form |7j: 

V = -V/ T ,, + 9r,h, 

h = -h/r h + g h gS, (1) 
S=(r- S)/t s - g s r)h. 

Here a point means a derivative in time r n , 77,, Tg — relaxation times of an order 
parameter rj(t), a conjugating field h(t) and a control parameter S(t), accord- 
ingly; g v , gh, gs - positive feed-back constants; r - pump intensity, measures the 
influence of environment. First elements in the left hand of the systems ([!} take 
into account dissipation effects, peculiar to the synergetic systems. Connection 
between the order parameter and conjugating field is linear (first equation), in 
that time as an evolution of the conjugating field h(t), and control parameter 
S(t) sets due to nonlinear feed-backs relations (second and third equations, re- 
spectively) . Principally that positive feed-backs which are provided by constant 
g v and gh result in an increase in conjugating field. These positive feed-backs 
are compensated by negative one due to principle of Le-Shatel'e. As a result 
one has decreasing in control parameter (see third equation in Eq. ([!])). 

Let us start the analysis of the system ([!} with passing to dimensionless 
variables. Such transition is arrived due to measuring of time t, order parameter 
rj, conjugating field h, and control parameter S in the followings units: 

iai"r, Ve oc (ghgs)~ 1/2 i 
h e oc (g^ghgs) 1/2 , S e oc {g n gh)~ x - 
In this time, dropping indexes, the system ([TJ becomes a form 

r) = — 77 + h, 

ah= -h + r)S, (2) 
eS — (r — S) — gh, 
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where a = Th/r^, e = ts/t v . The system ^ is written in supposition of linear 
dependence for order parameter relaxation time as 7^(77) = const. However, 
most real physical systems are characterized by the nonlinear relaxation pro- 
cesses. In this connection let us suppose that the order parameter relaxation 
time t, ; increasing with increase in order parameter r\ due to relation |10| : 

where k - positive constant which play a role of an dissipation intensity. From 
Eq. |3]) it is seen that relaxation time t v (if) is independent of order parameter 
sign. Except for that, relation ([3]) has practical application namely, it designs the 
action of optical filter, entered into the Fabry-Perot cavity of optically bistable 
system (for example solid-state laser). Such acting provides establishment of the 
stable periodic radiation (or time dissipative structure realization) | 11| . Using 
a dependence ^ , first equation of ^ is generalized by an additional nonlinear 
term f K = -(kt])/(1 + rf). 

If one consider system in external field, then one need to take into considera- 
tion external perturbations. In this article we shall model such perturbations by 
the external potential V e . Due to the standard catastrophe theory such potential 



is given by three types of catastrophes 12 . In general case one has 



V e = An + -if + -n 3 + -if + -n 5 , (4) 

where A, B, C, D, E - parameters of the theory. For a catastrophe A2 one has 
B = D = E = 0, for a catastrophe A3: C = E = and for a catastrophe A4: 
D = 0. The modified Lorenz-Haken system has the form 

f) = -V + h + f K (v) + feiv), 

h = -h + r/S, (5) 
S = (r — S) — rjh, 

where we suppose a ~ e ~ 1 and f e (j]) = — dV e /drj. Variation in parameters 
of / K (?7) and f e (r]) can to induce changing of the attractor topology in phase 
space. 



3 Chaos in a modified Lorenz-Haken system 

The system Eq.(|5| with nonlinear dependence of order parameter relaxation 
time versus order parameter in a form ^ (k ^ 0) but with absence of additional 
perturbations (V e = 0) was considered in [11] . It was shown that in such a 
case the semirestricted domain of system parameters (pump intensity r and 
dissipation intensity k) for dissipative structure realization is formed. It was set 
that in a case of linear dependence for order parameter relaxation time versus 
its values (k = 0) chaotic regime is not realize. In addition, it was found the 
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Figure 1: A Lypunov map of the modified Lorenz-Haken system ^ at n = 25.0 
and A = 0.1 



chaos domain, and it was determined conditions of chaos control. Finally it 
was defined fractal and statistical properties of corresponding chaotic strange 
attractor. 

The main goal in this work is to study an influence of external perturbation, 
on a regimes of transition to chaos in Lorenz-Haken system, generalized by non- 
linear relaxation time of order parameter in a form ([3]) . For external perturba- 
tions we shall use a potential of the fold catastrophe A 2 , i.e. V e = Arj+ l/ZCrf . 
To indicate a chaotic dynamic we shall use a method of Lyapunov exponents 
which is provided by a Benettin algorithm 13 . Due to this algorithm each of 



Lyaponov exponent (number is defined by dimension of corresponding phase 
space) determines a speed of convergence/divergence of any two initially nearby 
trajectories in a fixed direction in corresponding phase space, starting from 
points v(t) and v(t'). The divergence/convergence of such trajectories is given 
by the dependence 5v(t) — 5v(to)e AMt , where Am is a maximal (global) Lya- 
punov exponent, which is defined due to relation |13j 
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A M = A(5(t )) = lim - In 

1 — >OQ 1 



Sv(t) 



5v(t ) 



Here one takes into account an upper limit and \\v\\ is a norm; v = {rj, h, S}; T 
- full time. One can conclude that in a case Aj < 0, i ~ 1,2,3, and accordingly 
Am < all of phase trajectories will coincide to fixed point (stable node or 
stable focus). At Aj < 0, Aj < and A M = A k = 0, i ^ j ' ^ k i,j,k = 1,2,3 
phase trajectories will lie down on a stable limit circle (dissipative structure). 
If Ai < 0, Aj = and Am = A& > 0, i ^ j ^ k i, j, k = 1, 2, 3 a dynamics of the 
system is chaotic. A Lyapunov map of the modified Lorenz-Haken system (|5| 
at Hi = 25.0 and A = 0.1 is shown in Fig[T] Here by gradation of grey color the 
value of maximal (global) Lyapunov exponent is shown versus pump intensity r 
and parameter of an external potential C . White color determines the domains 
of stable system behaviour (phase space is characterized by a fixed point - 
stable node or a stable focus) . Grey color marks the domains of time dissipative 
structure existence (more dark domains correspond to the larger number of 



4 



oscillation periods). Domains of chaos are shown by a black colour. From Figjl] 
it is seen, that in a case of nonlinear term Crf absence in the external potential 
V e the existence of the chaotic mode requires large values of pump intensity. In 
case C 7^ : \C\ ~ 0.1 the chaotic mode exists at r ~ k. 

Let us analyze a picture of dynamical regimes reconstruction for the system 
|5]) with C = and A = 0.1 in detail. Lyapunov map and corresponding de- 
pendence of maximal Lyapunov exponent versus pump intensity at fixed value 
of dissipation intensity k with characteristic phase portraits are shown in Fig(2] 
Here due to earlier denoted scenario values of maximal Lyapunov exponent 
are presented by gradation of grey color versus pump intensity and dissipation 
intensity. It is necessary to note that dark curves (between a and b, b and c in a 
Lyapunov map) in the domain of dissipation structure existence (grey area) are 
determine the parameters values of doubling period bifurcation. Below the map 
there is a dependence of maximal Lyapunov exponent versus pump intensity r 
at A = 0.1, C = 0.0 and k = 25.0. A presence of two pronounced peaks in 
the domain of zero values of maximal Lyapunov exponent (fluctuations around 
zero are connected with an error of numeral solution of systems (|5|) determines 
points of doubling period bifurcation. Corresponding phase portraits illustrat- 
ing dissipative structures with one, n and m > n periods are shown with the 
help of insertions a) , b) and c) , accordingly. It is principally, that as transition 
from a) to b) , as transition from b) to c) characterizes by appearance of a few 
additional harmonics. As it is seen from dependence Ajw(r) in a point d maxi- 
mal Lyapunov exponent has positive value and phase space characterizes by the 
irregular behaviour of trajectory. Increasing in pump intensity leads to dissipa- 
tive structure formation with one period (cf. phase portraits d) and e)). Next 
increasing in r leads to positive values of Am and phase space characterizes by 
chaos existence (corresponding phase portrait is shown with the help of insertion 
f). So, one can conclude that at A — 0.1, C — 0.0 and k = 25.0 with an increase 
in pump intensity r a transition to chaotic regime occurs due to Ruelle-Takens 
scenario, when only negligible number of doubling period bifurcations leads to 
chaos 1 14] , With a decreasing in r (from r > r c to r < r c ) maximal Lyapunov 
exponent takes positive value at r = r c in spontaneous manner. Corresponding 
transition to chaos takes a place through intermittency flB"] . 

Next, let us consider the domain of chaos, shown in Figjljat C ^ 0. Lya- 
punov map at A = 0.1 and C = 0.1 is shown in Fig|3] Below the map as well 
as in previous case a dependence of maximal Lyapunov exponent versus pump 
intensity at A = 0.1, C = 0.1 and n = 25.0 and corresponding phase portraits 
are shown. Unlike to the previous case here with an increase in pump intensity 
r a successive complication of attractor due to doubling period bifurcation is 
observed (see corresponding phase portraits in points a, b, c, d and e). Thus, 
in such a case (A = 0.1, C = 0.1 and n = 25.0) an increasing in r results to 
transition to chaos due to Feigenbaum scenario [M] . Chaotic attractor is shown 
with the help of insertion f). As in a previous case a decrease in r leads to 
transition to chaos at r = r c through intermittency |15| . 

It is well known that dynamical systems can realize four types of attractors 
in phase space, namely: non chaotic non strange attractor, chaotic non strange 
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Figure 2: Lyapunov map of the modified Lorenz-Haken system ([5| at A = 0.1 
and C — 0.0; dependence of maximal Lyapunov exponent at A — 0.1, C = 0.0 
and k = 25.0 and corresponding phase portraits 



attractor, strange non chaotic attractor and chaotic strange attractor. So, for 
chaotic attractor one has (Am > 0) and for an strange one a fractal dimension 
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Figure 3: Lyapunov map of the modified Lorenz-Haken system ([5| at A = 0.1 
and C — 0.1; dependence of maximal Lyapunov exponent at A — 0.1, C = 0.1 
and k = 25.0 and corresponding phase portraits 



of an attractor is fractional. In [16] it was shown that the fractal dimension of 
an attractor which is realized in the dynamical system, is determined with the 
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help of Lyapunov exponents due to relation 



D~ 1.5 + 0.5W1-8 — — , 

V ^min 

where A m j„ is a minimal Lyapunov exponent. Thus, for considered attractor 
in a point f at k = 25.0, r = 98.0, A = 0.1 and C = (see Figj2j one has: 
A M = 0.2774, A rnin = -16.948 and, accordingly, D ~ 2.032. For the attractor 
at k = 25.0, r = 34.5, A = 0.1 and C = 0.1 (see Fig(3} one has A M = 0.16545, 
A m m = —3.631 and accordingly, D ~ 2.084. Thus, attractors in Figj2j" and 
Fig(3f are strange and chaotic. 



4 Conclusions 

We have studied an influence of nonlinear dissipation and external perturba- 
tions onto transition scenarious to chaos in Lorenz-Haken system. Dissipation 
processes are defined due to nonlinear dependence of order parameter relaxation 
time r n versus its values. External perturbations are modeled by a potential of 
the fold catastrophe A2. 

From a physical view point we have considered an absorptive optical bista- 
bility system 17] . At that time used nonlinear dissipation relation related to the 
possibility of the additional medium in the Fabry-Perot cavity (phthalocyanine 
fluid, gases SFq, BaCl^, and CO2 (7|) to absorbing signals with weak intensi- 
ties. Meanwhile, external perturbations model an influence of optical modulator 
which sets additional interphotons interaction processes in the Fabry-Perot cav- 
ity. 

We have considered a case of commcnsurability of relaxation times for order 
parameter, conjugated field and control parameter. We have shown that varying 
in external potential parameters values leads to parameters domain formation 
with r ~ k and with r by order of magnitude greater than k of chaos realization. 
In considered system transitions from regular to chaotic dynamics can be of 
Ruelle-Takens scenario, Feigenbaum scenario, or through intermittency. 
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